Sex differences in DNA methylation across gestation: a large scale, cross-cohort, multi-tissue analysis.

Biological sex is a key variable influencing many physiological systems. Disease prevalence as well as treatment success can be modified by sex. Differences emerge already early in life and include pregnancy complications and adverse birth outcomes. The placenta is a critical organ for fetal development and shows sex-based differences in the expression of hormones and cytokines. Epigenetic regulation, such as DNA methylation (DNAm), may underlie the previously reported placental sexual dimorphism. We associated placental DNAm with fetal sex in three cohorts. Individual cohort results were meta-analyzed with random-effects modelling. CpG-sites differentially methylated with sex were further investigated regarding pathway enrichment, overlap with methylation quantitative trait loci (meQTLs), and hits from phenome-wide association studies (PheWAS). We evaluated the consistency of findings across tissues (CVS, i.e. chorionic villus sampling from early placenta, and cord blood) as well as with gene expression. We identified 10,320 epigenome-wide significant sex-differentially methylated probes (DMPs) spread throughout the epigenome of the placenta at birth. Most DMPs presented with lower DNAm levels in females. DMPs mapped to genes upregulated in brain, were enriched for neurodevelopmental pathways and significantly overlapped with meQTLs and PheWAS hits. Effect sizes were moderately correlated between CVS and placenta at birth, but only weakly correlated between birth placenta and cord blood. Sex differential gene expression in birth placenta was less pronounced and implicated genetic regions only marginally overlapped with those associated with differential DNAm. Our study provides an integrative perspective on sex-differential DNAm in perinatal tissues underscoring the possible link between placenta and brain.

Genetic predisposition for negative affect predicts mental health burden during the COVID-19 pandemic.

The coronavirus disease 2019 (COVID-19) pandemic was accompanied by an increase in mental health challenges including depression, stress, loneliness, and anxiety. Common genetic variants can contribute to the risk for psychiatric disorders and may present a risk factor in times of crises. However, it is unclear to what extent polygenic risk played a role in the mental health response to the COVID-19 pandemic. In this study, we investigate whether polygenic scores (PGSs) for mental health-related traits can distinguish between four resilience-vulnerability trajectories identified during the COVID-19 pandemic and associated lockdowns in 2020/21. We used multinomial regression in a genotyped subsample (n = 1316) of the CovSocial project. The most resilient trajectory characterized by the lowest mental health burden and the highest recovery rates served as the reference group. Compared to this most resilient trajectory, a higher value on the PGS for the well-being spectrum decreased the odds for individuals to be in one of the more vulnerable trajectories (adjusted R-square = 0.3%). Conversely, a higher value on the PGS for neuroticism increased the odds for individuals to be in one of the more vulnerable trajectories (adjusted R-square = 0.2%). Latent change in mental health burden extracted from the resilience-vulnerability trajectories was not associated with any PGS. Although our findings support an influence of PGS on mental health during COVID-19, the small added explained variance suggests limited utility of such genetic markers for the identification of vulnerable individuals in the general population.

Serial Diurnal Salivary Cortisol Profiles in 667 pregnant women - Association with Cardiometabolic Complications.

CONTEXT: Maternal obesity, hypertensive pregnancy disorders and gestational diabetes (GDM) are linked to an increased risk of negative offspring health outcomes. This association may be mediated by maternal hypothalamic-pituitary-adrenal axis (HPA axis) activity, resulting in elevated maternal cortisol levels and fetal exposure, but evidence remains scarce. OBJECTIVE: We examined (1) maternal diurnal cortisol profiles longitudinally across gestation, and (2) explored associations with maternal cardiometabolic complications. DESIGN: Women in the InTraUterine sampling in early pregnancy (ITU) study (n=667) provided seven salivary cortisol samples from awakening to bedtime up to three times during pregnancy (median gestational week 19.3, 25.7, and 38.1, n=9,356 samples). Changes in cortisol awakening response and diurnal slope (indicative of HPA-axis activity) and their associations with maternal body mass index (BMI), hypertensive pregnancy disorders and GDM were examined using linear mixed models. RESULTS: The cortisol awakening response declined in in 60%-67% of women, and the diurnal slope attenuated from early to late pregnancy (b = 0.006, p = .001). Higher BMI was associated with less decline in cortisol awakening response (b= 0.031, p = .0004), and less attenuation in diurnal slope from early to late pregnancy (b = -0.001, p = .006). Hypertensive pregnancy disorders and GDM were not significantly associated with diurnal cortisol profiles. CONCLUSIONS: The attenuation in cortisol awakening response and diurnal slope support HPA-axis hypo-responsivity during pregnancy. Less attenuation of both markers in women with a higher BMI may indicate reduced adaption of the HPA-axis to pregnancy, presenting a mechanistic link to offspring health outcomes.

Cohort profile: BioMD-Y (biopsychosocial factors of major depression in youth) - a biobank study on the molecular genetics and environmental factors of depression in children and adolescents in Munich.

PURPOSE: BioMD-Y is a comprehensive biobank study of children and adolescents with major depression (MD) and their healthy peers in Germany, collecting a host of both biological and psychosocial information from the participants and their parents with the aim of exploring genetic and environmental risk and protective factors for MD in children and adolescents. PARTICIPANTS: Children and adolescents aged 8-18 years are recruited to either the clinical case group (MD, diagnosis of MD disorder) or the typically developing control group (absence of any psychiatric condition). FINDINGS TO DATE: To date, four publications on both genetic and environmental risk and resilience factors (including FKBP5, glucocorticoid receptor activation, polygenic risk scores, psychosocial and sociodemographic risk and resilience factors) have been published based on the BioMD-Y sample. FUTURE PLANS: Data collection is currently scheduled to continue into 2026. Research questions will be further addressed using available measures.

Human cortical neurogenesis is altered via glucocorticoid-mediated regulation of ZBTB16 expression.

Glucocorticoids are important for proper organ maturation, and their levels are tightly regulated during development. Here, we use human cerebral organoids and mice to study the cell-type-specific effects of glucocorticoids on neurogenesis. We show that glucocorticoids increase a specific type of basal progenitors (co-expressing PAX6 and EOMES) that has been shown to contribute to cortical expansion in gyrified species. This effect is mediated via the transcription factor ZBTB16 and leads to increased production of neurons. A phenome-wide Mendelian randomization analysis of an enhancer variant that moderates glucocorticoid-induced ZBTB16 levels reveals causal relationships with higher educational attainment and altered brain structure. The relationship with postnatal cognition is also supported by data from a prospective pregnancy cohort study. This work provides a cellular and molecular pathway for the effects of glucocorticoids on human neurogenesis that relates to lasting postnatal phenotypes.

Identifying genetic differences between bipolar disorder and major depression through multiple GWAS.

Background: Accurate diagnosis of bipolar disorder (BD) is difficult in clinical practice, with an average delay between symptom onset and diagnosis of about 7 years. A key reason is that the first manic episode is often preceded by a depressive one, making it difficult to distinguish BD from unipolar major depressive disorder (MDD). Aims: Here, we use genome-wide association analyses (GWAS) to identify differential genetic factors and to develop predictors based on polygenic risk scores that may aid early differential diagnosis. Methods: Based on individual genotypes from case-control cohorts of BD and MDD shared through the Psychiatric Genomics Consortium, we compile case-case-control cohorts, applying a careful merging and quality control procedure. In a resulting cohort of 51,149 individuals (15,532 BD cases, 12,920 MDD cases and 22,697 controls), we perform a variety of GWAS and polygenic risk scores (PRS) analyses. Results: While our GWAS is not well-powered to identify genome-wide significant loci, we find significant SNP-heritability and demonstrate the ability of the resulting PRS to distinguish BD from MDD, including BD cases with depressive onset. We replicate our PRS findings, but not signals of individual loci in an independent Danish cohort (iPSYCH 2015 case-cohort study, N=25,966). We observe strong genetic correlation between our case-case GWAS and that of case-control BD. Conclusions: We find that MDD and BD, including BD with a depressive onset, are genetically distinct. Further, our findings support the hypothesis that Controls - MDD - BD primarily lie on a continuum of genetic risk. Future studies with larger and richer samples will likely yield a better understanding of these findings and enable the development of better genetic predictors distinguishing BD and, importantly, BD with depressive onset from MDD.

Epigenetics of prenatal stress in humans: the current research landscape.

Fetal exposure to prenatal stress can have significant consequences on short- and long-term health. Epigenetic mechanisms, especially DNA methylation (DNAm), are a possible process how these adverse environmental events could be biologically embedded. We evaluated candidate gene as well as epigenome-wide association studies associating prenatal stress and DNAm changes in peripheral tissues; however, most of these findings lack robust replication. Prenatal stress-associated epigenetic changes have also been linked to child health including internalizing problems, neurobehavioral outcomes and stress reactivity. Future studies should focus on refined measurement and definition of prenatal stress and its timing, ideally also incorporating genomic as well as longitudinal information. This will provide further opportunities to enhance our understanding of the biological embedding of prenatal stress exposure.

DNA methylation patterns of FKBP5 regulatory regions in brain and blood of humanized mice and humans.

Humanized mouse models can be used to explore human gene regulatory elements (REs), which frequently lie in non-coding and less conserved genomic regions. Epigenetic modifications of gene REs, also in the context of gene x environment interactions, have not yet been explored in humanized mouse models. We applied high-accuracy measurement of DNA methylation (DNAm) via targeted bisulfite sequencing (HAM-TBS) to investigate DNAm in three tissues/brain regions (blood, prefrontal cortex and hippocampus) of mice carrying the human FK506-binding protein 5 (FKBP5) gene, an important candidate gene associated with stress-related psychiatric disorders. We explored DNAm in three functional intronic glucocorticoid-responsive elements (at introns 2, 5, and 7) of FKBP5 at baseline, in cases of differing genotype (rs1360780 single nucleotide polymorphism), and following application of the synthetic glucocorticoid dexamethasone. We compared DNAm patterns in the humanized mouse (N = 58) to those in human peripheral blood (N = 447 and N = 89) and human postmortem brain prefrontal cortex (N = 86). Overall, DNAm patterns in the humanized mouse model seem to recapitulate DNAm patterns observed in human tissue. At baseline, this was to a higher extent in brain tissue. The animal model also recapitulated effects of dexamethasone on DNAm, especially in peripheral blood and to a lesser extent effects of genotype on DNAm. The humanized mouse model could thus assist in reverse translation of human findings in psychiatry that involve genetic and epigenetic regulation in non-coding elements.

A Pregnancy and Childhood Epigenetics Consortium (PACE) meta-analysis highlights potential relationships between birth order and neonatal blood DNA methylation.

Higher birth order is associated with altered risk of many disease states. Changes in placentation and exposures to in utero growth factors with successive pregnancies may impact later life disease risk via persistent DNA methylation alterations. We investigated birth order with Illumina DNA methylation array data in each of 16 birth cohorts (8164 newborns) with European, African, and Latino ancestries from the Pregnancy and Childhood Epigenetics Consortium. Meta-analyzed data demonstrated systematic DNA methylation variation in 341 CpGs (FDR adjusted P < 0.05) and 1107 regions. Forty CpGs were located within known quantitative trait loci for gene expression traits in blood, and trait enrichment analysis suggested a strong association with immune-related, transcriptional control, and blood pressure regulation phenotypes. Decreasing fertility rates worldwide with the concomitant increased proportion of first-born children highlights a potential reflection of birth order-related epigenomic states on changing disease incidence trends.

Maternal educational attainment in pregnancy and epigenome-wide DNA methylation changes in the offspring from birth until adolescence.

Maternal educational attainment (MEA) shapes offspring health through multiple potential pathways. Differential DNA methylation may provide a mechanistic understanding of these long-term associations. We aimed to quantify the associations of MEA with offspring DNA methylation levels at birth, in childhood and in adolescence. Using 37 studies from high-income countries, we performed meta-analysis of epigenome-wide association studies (EWAS) to quantify the associations of completed years of MEA at the time of pregnancy with offspring DNA methylation levels at birth (n = 9 881), in childhood (n = 2 017), and adolescence (n = 2 740), adjusting for relevant covariates. MEA was found to be associated with DNA methylation at 473 cytosine-phosphate-guanine sites at birth, one in childhood, and four in adolescence. We observed enrichment for findings from previous EWAS on maternal folate, vitamin-B12 concentrations, maternal smoking, and pre-pregnancy BMI. The associations were directionally consistent with MEA being inversely associated with behaviours including smoking and BMI. Our findings form a bridge between socio-economic factors and biology and highlight potential pathways underlying effects of maternal education. The results broaden our understanding of bio-social associations linked to differential DNA methylation in multiple early stages of life. The data generated also offers an important resource to help a more precise understanding of the social determinants of health.

High-throughput screening of glucocorticoid-induced enhancer activity reveals mechanisms of stress-related psychiatric disorders.

Exposure to stressful life events increases the risk for psychiatric disorders. Mechanistic insight into the genetic factors moderating the impact of stress can increase our understanding of disease processes. Here, we test 3,662 single nucleotide polymorphisms (SNPs) from preselected expression quantitative trait loci in massively parallel reporter assays to identify genetic variants that modulate the activity of regulatory elements sensitive to glucocorticoids, important mediators of the stress response. Of the tested SNP sequences, 547 were located in glucocorticoid-responsive regulatory elements of which 233 showed allele-dependent activity. Transcripts regulated by these functional variants were enriched for those differentially expressed in psychiatric disorders in the postmortem brain. Phenome-wide Mendelian randomization analysis in 4,439 phenotypes revealed potentially causal associations specifically in neurobehavioral traits, including major depression and other psychiatric disorders. Finally, a functional gene score derived from these variants was significantly associated with differences in the physiological stress response, suggesting that these variants may alter disease risk by moderating the individual set point of the stress response.

Transdiagnostic evaluation of epigenetic age acceleration and burden of psychiatric disorders.

Different psychiatric disorders as well as exposure to adverse life events have individually been associated with multiple age-related diseases and mortality. Age acceleration in different epigenetic clocks can serve as biomarker for such risk and could help to disentangle the interplay of psychiatric comorbidity and early adversity on age-related diseases and mortality. We evaluated five epigenetic clocks (Horvath, Hannum, PhenoAge, GrimAge and DunedinPoAm) in a transdiagnostic psychiatric sample using epigenome-wide DNA methylation data from peripheral blood of 429 subjects from two studies at the Max Planck Institute of Psychiatry. Burden of psychiatric disease, represented by a weighted score, was significantly associated with biological age acceleration as measured by GrimAge and DunedinPoAm (R2-adj. 0.22 and 0.33 for GrimAge and DunedinPoAm, respectively), but not the other investigated clocks. The relation of burden of psychiatric disease appeared independent of differences in socioeconomic status and medication. Our findings indicate that increased burden of psychiatric disease may associate with accelerated biological aging. This highlights the importance of medical management of patients with multiple psychiatric comorbidities and the potential usefulness of specific epigenetic clocks for early detection of risk and targeted intervention to reduce mortality in psychiatric patients.

Epigenome-wide meta-analysis of prenatal maternal stressful life events and newborn DNA methylation.

Prenatal maternal stressful life events are associated with adverse neurodevelopmental outcomes in offspring. Biological mechanisms underlying these associations are largely unknown, but DNA methylation likely plays a role. This meta-analysis included twelve non-overlapping cohorts from ten independent longitudinal studies (N = 5,496) within the international Pregnancy and Childhood Epigenetics consortium to examine maternal stressful life events during pregnancy and DNA methylation in cord blood. Children whose mothers reported higher levels of cumulative maternal stressful life events during pregnancy exhibited differential methylation of cg26579032 in ALKBH3. Stressor-specific domains of conflict with family/friends, abuse (physical, sexual, and emotional), and death of a close friend/relative were also associated with differential methylation of CpGs in APTX, MyD88, and both UHRF1 and SDCCAG8, respectively; these genes are implicated in neurodegeneration, immune and cellular functions, regulation of global methylation levels, metabolism, and schizophrenia risk. Thus, differences in DNA methylation at these loci may provide novel insights into potential mechanisms of neurodevelopment in offspring.

Associations of psychiatric disease and ageing with FKBP5 expression converge on superficial layer neurons of the neocortex.

Identification and characterisation of novel targets for treatment is a priority in the field of psychiatry. FKBP5 is a gene with decades of evidence suggesting its pathogenic role in a subset of psychiatric patients, with potential to be leveraged as a therapeutic target for these individuals. While it is widely reported that FKBP5/FKBP51 mRNA/protein (FKBP5/1) expression is impacted by psychiatric disease state, risk genotype and age, it is not known in which cell types and sub-anatomical areas of the human brain this occurs. This knowledge is critical to propel FKBP5/1-targeted treatment development. Here, we performed an extensive, large-scale postmortem study (n = 1024) of FKBP5/1, examining neocortical areas (BA9, BA11 and ventral BA24/BA24a) derived from subjects that lived with schizophrenia, major depression or bipolar disorder. With an extensive battery of RNA (bulk RNA sequencing, single-nucleus RNA sequencing, microarray, qPCR, RNAscope) and protein (immunoblot, immunohistochemistry) analysis approaches, we thoroughly investigated the effects of disease state, ageing and genotype on cortical FKBP5/1 expression including in a cell type-specific manner. We identified consistently heightened FKBP5/1 levels in psychopathology and with age, but not genotype, with these effects strongest in schizophrenia. Using single-nucleus RNA sequencing (snRNAseq; BA9 and BA11) and targeted histology (BA9, BA24a), we established that these disease and ageing effects on FKBP5/1 expression were most pronounced in excitatory superficial layer neurons of the neocortex, and this effect appeared to be consistent in both the granular and agranular areas examined. We then found that this increase in FKBP5 levels may impact on synaptic plasticity, as FKBP5 gex levels strongly and inversely correlated with dendritic mushroom spine density and brain-derived neurotrophic factor (BDNF) levels in superficial layer neurons in BA11. These findings pinpoint a novel cellular and molecular mechanism that has potential to open a new avenue of FKBP51 drug development to treat cognitive symptoms in psychiatric disorders.

Maternal exposure to childhood maltreatment and mental and behavioral disorders in children.

Exposure to maltreatment in childhood is associated with lifelong risk of mental and behavioral disorders. Whether the effects extend to the next generation remains unclear. We examined whether maternal exposure to childhood abuse and neglect in her own childhood were associated with mental and behavioral disorders and psychiatric symptoms in her children, and whether maternal lifetime mental and behavioral disorders or lower education level mediated or added to the effects. Mothers (n = 2252) of the Prediction and Prevention of Preeclampsia and Intrauterine Growth Restriction cohort study completed the Childhood Trauma Questionnaire and reported on their education and their 7.0-12.1-year-old children's psychiatric symptoms using the Strengths and Difficulties Questionnaire. We identified lifetime mental and behavioral disorder diagnoses for the mothers and diagnoses for their children from birth (2006-2010) until 8.4-12.8 years (12/31/2018) from Care Register for Health Care. We found that maternal exposure to childhood abuse, but not neglect, was associated with higher hazards of mental and behavioral disorders (hazard ratio 1.20, 95% confidence interval 1.06-1.37) in children. These associations were partially mediated by maternal mental and behavioral disorders and education (proportion of effect size mediated: 23.8% and 15.1%, respectively), which together with maternal exposure to childhood abuse added to the hazard of mental and behavioral disorders in children. Similar associations were found for maternal exposure to childhood abuse and neglect with psychiatric symptoms in children. To conclude, maternal exposure to childhood maltreatment is associated with mental and behavioral disorders and psychiatric symptoms in children. Our findings call for interventions to prevent intergenerational transmission.

A meta-analysis of pre-pregnancy maternal body mass index and placental DNA methylation identifies 27 CpG sites with implications for mother-child health.

Higher maternal pre-pregnancy body mass index (ppBMI) is associated with increased neonatal morbidity, as well as with pregnancy complications and metabolic outcomes in offspring later in life. The placenta is a key organ in fetal development and has been proposed to act as a mediator between the mother and different health outcomes in children. The overall aim of the present work is to investigate the association of ppBMI with epigenome-wide placental DNA methylation (DNAm) in 10 studies from the PACE consortium, amounting to 2631 mother-child pairs. We identify 27 CpG sites at which we observe placental DNAm variations of up to 2.0% per 10 ppBMI-unit. The CpGs that are differentially methylated in placenta do not overlap with CpGs identified in previous studies in cord blood DNAm related to ppBMI. Many of the identified CpGs are located in open sea regions, are often close to obesity-related genes such as GPX1 and LGR4 and altogether, are enriched in cancer and oxidative stress pathways. Our findings suggest that placental DNAm could be one of the mechanisms by which maternal obesity is associated with metabolic health outcomes in newborns and children, although further studies will be needed in order to corroborate these findings.

Genome-Wide Copy Number Variant and High-Throughput Transcriptomics Analyses of Placental Tissues Underscore Persisting Child Susceptibility in At-Risk Pregnancies Cleared in Standard Genetic Testing.

Several studies have shown that children from pregnancies with estimated first-trimester risk based on fetal nuchal translucency thickness and abnormal maternal serum pregnancy protein and hormone levels maintain a higher likelihood of adverse outcomes, even if initial testing for known genetic conditions is negative. We used the Finnish InTraUterine cohort (ITU), which is a comprehensively characterized perinatal cohort consisting of 943 mothers and their babies followed throughout pregnancy and 18 months postnatally, including mothers shortlisted for prenatal genetic testing but cleared for major aneuploidies (cases: n = 544, 57.7%) and control pregnancies (n = 399, 42.3%). Using genome-wide genotyping and RNA sequencing of first-trimester and term placental tissue, combined with medical information from registry data and maternal self-report data, we investigated potential negative medical outcomes and genetic susceptibility to disease and their correlates in placenta gene expression. Case mothers did not present with higher levels of depression, perceived stress, or anxiety during pregnancy. Case children were significantly diagnosed more often with congenital malformations of the circulatory system (4.12 (95% CI [1.22−13.93]) higher hazard) and presented with significantly more copy number duplications as compared to controls (burden analysis, based on all copy number variants (CNVs) with at most 10% frequency, 823 called duplications in 297 cases versus 626 called duplications in 277 controls, p = 0.01). Fifteen genes showed differential gene expression (FDR < 0.1) in association with congenital malformations in first-trimester but not term placenta. These were significantly enriched for genes associated with placental dysfunction. In spite of normal routine follow-up prenatal testing results in early pregnancy, case children presented with an increased likelihood of negative outcomes, which should prompt vigilance in follow-up during pregnancy and after birth.

Effects of stressful life-events on DNA methylation in panic disorder and major depressive disorder.

BACKGROUND: Panic disorder (PD) is characterized by recurrent panic attacks and higher affection of women as compared to men. The lifetime prevalence of PD is about 2-3% in the general population leading to tremendous distress and disability. Etiologically, genetic and environmental factors, such as stress, contribute to the onset and relapse of PD. In the present study, we investigated epigenome-wide DNA methylation (DNAm) in respond to a cumulative, stress-weighted life events score (wLE) in patients with PD and its boundary to major depressive disorder (MDD), frequently co-occurring with symptoms of PD. METHODS: DNAm was assessed by the Illumina HumanMethylation450 BeadChip. In a meta-analytic approach, epigenome-wide DNAm changes in association with wLE were first analyzed in two PD cohorts (with a total sample size of 183 PD patients and 85 healthy controls) and lastly in 102 patients with MDD to identify possible overlapping and opposing effects of wLE on DNAm. Additionally, analysis of differentially methylated regions (DMRs) was conducted to identify regional clusters of association. RESULTS: Two CpG-sites presented with p-values below 1 × 10-05 in PD: cg09738429 (p = 6.40 × 10-06, located in an intergenic shore region in next proximity of PYROXD1) and cg03341655 (p = 8.14 × 10-06, located in the exonic region of GFOD2). The association of DNAm at cg03341655 and wLE could be replicated in the independent MDD case sample indicating a diagnosis independent effect. Genes mapping to the top hits were significantly upregulated in brain and top hits have been implicated in the metabolic system. Additionally, two significant DMRs were identified for PD only on chromosome 10 and 18, including CpG-sites which have been reported to be associated with anxiety and other psychiatric phenotypes. CONCLUSION: This first DNAm analysis in PD reveals first evidence of small but significant DNAm changes in PD in association with cumulative stress-weighted life events. Most of the top associated CpG-sites are located in genes implicated in metabolic processes supporting the hypothesis that environmental stress contributes to health damaging changes by affecting a broad spectrum of systems in the body.

Early adversity as the prototype gene × environment interaction in mental disorders?

Childhood adversity (CA) as a significant stressor has consistently been associated with the development of mental disorders. The interaction between CA and genetic variants has been proposed to play a substantial role in disease etiology. In this review, we focus on the gene by environment (GxE) paradigm, its background and interpretation and stress the necessity of its implementation in psychiatric research. Further, we discuss the findings supporting GxCA interactions, ranging from candidate gene studies to polygenic and genome-wide approaches, their strengths and limitations. To illustrate potential underlying epigenetic mechanisms by which GxE effects are translated, we focus on results from FKBP5 × CA studies and discuss how molecular evidence can supplement previous GxE findings. In conclusion, while GxE studies constitute a valuable line of investigation, more harmonized GxE studies in large, deep-phenotyped, longitudinal cohorts, and across different developmental stages are necessary to further substantiate and understand reported GxE findings.